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STEADY  STATE  MAGNETIC  DIFFUSION  FROM  RESISTIVE  INTERCHANGE  MODES 

It  is  now  generally  assumed  that  reverse  field  pinches  (RFPs) 

can  exist  neither  in  states  which  are  MHD  unstable  nor  in  states  which 

are  tearing  mode  unstable.  Indeed,  diffuse  pinch  profiles  which  are 

12  3 

stable  to  all  of  these  have  been  calculated.  *  *  It  seems  reasonable 

that  this  conclusion  also  applies  to  spheromaks.  The  remaining 

problem  is  the  effect  of  resistive  interchange  modes.  The  simplest 

theory,  which  mocks  up  pressure  gradient  with  an  effective  gravity, 

4  5  6 

indicates  that  any  pressure  gradient  drives  these  modes  unstable.  *  * 

Apparently  the  plasma  must  either  exist  in  such  an  unstable  state,  or 

else  be  driven  to  the  Taylor  configuration  with  zero  pressure 

gradient. ^  This  seems  to  be  a  crucial  issue  and  one  can  examine  it 

in  two  ways.  First  one  can  do  the  linear  theory  more  accurately, 

modeling  geometry  and/or  kinetic  effects  more  realistically,  and  hope 

that  stable  regimes  emerge.  Second,  one  can  examine  whether  the 

plasma  can  exist  in  the  presence  of  these  modes.  This  report  follows 

the  second  approach.  It  finds  that  the  plasma  can  indeed  exist  (but 

has  anomalous  transport),  with  a  spectrum  of  nonlinearly  saturated 

resistive  interchange  modes.  The  anomalous  magnetic  diffusion  is 

calculated  in  terms  of  the  fluid  velocity  fluctuation  and  an  estimate 

is  then  given  of  an  upperbound  for  this.  Also,  it  is  worth  pointing 

out  that  since  mode  rational  surfaces  for  resistive  g  modes  occur 
Note:  Manuscript  submitted  January  17,  1980. 
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everywhere  in  an  RFP,  a  theory  which  predicts  stabilization  by 
flattening  the  pressure  profile  over  for  instance  a  few  widths  of  the 
interaction  region  is  not  a  viable  theory,  because  these  interaction 
regions  will  almost  certainly  overlap. 

This  work  is  motivated  to  a  large  extent  by  recent  studies  in 

Q 

Eta-Beta  II.  There  it  was  found  that  a  reverse  field  pinch  plasma 
14 

with  n  a  2  x  10  ,  T  «  100  ev  and  I  =  2Q0KA  could  exist  in  a  quiescent 
state  for  as  long  as  500  y  sec  before  disrupting.  The  growth  time  for 
the  resistive  g  mode  is  less  than  10  y  sec,  while  the  resistive  diffusion 
time  is  about  10-20  milisec.  Thus  something  allows  this  plasma  to 
exist  for  many  growth  times,  but  to  be  lost  rapidly  compared  to  a 
classical  diffusion  time.  Also,  while  magnetic  probes  show  that  low 
frequency  fluctuations  virtually  disappear  during  the  quiescent  phase, 
the  high  frequency  fluctuations  are  reduced,  but  are  still  present. 

During  this  quiescent  period,  the  electron  temperature  increases  by 
50-100  ev.  It  is  simple  to  show  that  classical  Ohmic  heating  will 
give  rise  to  a  much  larger  temperature  increase,  and  also  to  show  that 
ion  thermal  conduction  does  not  significantly  cool  the  plasma  during 
a  relevant  time  scale.  Thus  there  appears  to  be  an  anomalous  heat  loss 
also. 

There  are  several  reasons  to  examine  the  resistive  inter¬ 
change  node  as  well  as  other  micro-instabilities.  First,  the  high 
frequency  probes  pick  up  magnetic  signals  so  the  mode  cannot  be 


2 


t 


purely  electrostatic.  Secondly,  an  RFP  has  strong  shear  and  this 
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usually  has  a  great  stabilizing  effect  on  micro-instabilities, 
whereas  a  resistive  interchange  mode  cannot  be  shear  stabilized. 

Finally  a  resistive  interchange  mode  is  simpler  than  a  micro-instability 
since  it  can  be  described  within  a  fluid  framework,  so  a  nonlinear 
theory  of  it  is  interesting  in  its  own  right. 

We  begin  by  briefly  reviewing  (following  Ref.  6)  the  linear 
thes>ry  of  the  resistive  interchange  mode.  The  idea  is  that  resistivity 
allows  the  fluid  to  slip  through  the  field  lines,  and  thereby  defeat 
the  shear  stabilization  which  would  occur  in  ideal  MHD.  That  is,  the 
V  x  B  force  in  the  plasma  does  not  generate  an  electric  field,  but  is 
balanced  by  resistive  heating,  so  that  instead  of  E  +  ^  x  B  ■=  0, 

Ohm's  law  becomes 

n  3  *  J  x  B  (1) 

where  a  top  squiggle  indicates  a  fluctuating  quantity,  and  n  is  assumed 
to  be  a  constant.  In  slab  geometry,  with  B=B(i  +  r~  i  ),  incom- 

O  -Z  L  “Y 
S 


pressible  perturbed  motion  and  with  the  fluctuating  quanties  varying  as 
f(x)  exp  (iky  +  yt)  +  c.c..  Ref.  6  calculates  the  mode  structure  and 
growth  rate.  The  result  is 


T 

where  g  is  the  equivalent  gravity  g  'v  xL,  and  R  is  the  radius  of 

MK 

curvature  of  the  field  line.  The  plasma  is  unstable  only  if  and 

R  have  opposite  sign.  For  a  reversed  field  pinch,  R  is  roughly  the 

minor  radius  r.  Also,  it  is  Important  to  note  that  for  a  resistive 

interchange  mode,  E  -  0  according  to  Eq.  (1).  Notice  that  all 

quantities  except  are  localized  to  within  a  distance  *.  of  the 

rational  surface.  Far  from  the  rational  surface,  is  actually 
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determined  by  the  full  MHD  equations  and  Eq.  (2d)  is  not  accurate, 
however,  this  is  not  important  in  what  follows.  In  a  future  publi¬ 
cation,  a  more  detailed  discussion  of  the  linear  theory  (and  nonlinear 
theory)  will  be  given  in  both  slab  and  more  realistic  geometry. 

We  now  turn  to  the  quasi-linear  theory.  In  ideal  MHD,  quasi- 
linear  theory  shows  that  the  background  plasma  responds  to  the 

unstable  mode  only  while  that  mode  is  growing. ^  This  corresponds  to 
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nonresonant  quasi-linear  theory  in  an  infinite  homogeneous  plasma. 

However,  for  a  resistive  interchange  mode,  dissipation  is  present. 

As  we  will  see,  this  means  that  the  background  plasma  can  respond  even 

if  the  fluctuations  are  at  steady  state.  This  corresponds  to  resonant 
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quasi-linear  theory  in  an  infinite  homogeneous  plasma. 
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Now  imagine  that  some  nonlinear  effect  stops  the  growth  of  the 
fluctuation  (i.e.,  so  y  *  0)  at  some  saturated  value  which  we  will 
denote  VQ.  (In  addition  to  specifying  a  Vq,  there  are  other  subtle 
requirements  on  such  a  nonlinear  effect  which  will  be  discussed  more 
fully  else  where) .  This  paper  does  not  speculate  on  what  this  non¬ 
linear  effect  is;  it  only  assumes  a  saturated  value  for  Vq  and  proceeds. 

The  idea  then  is  that  the  plasma  is  not  quiescent,  but  has  a 
fluctuating  velocity.  This  fluctuating  velocity  is  balanced  by  Ohmic 
dissipation  in  steady  state  according  to  Eq.  (1).  However,  this 
fluctuating  velocity  gives  rise  to  a  steady  state  electric  field  in 
the  z  direction  according  to  Ohms  law: 

E  «  -  -  fv  B  *-  V  B  *^+C*C*  +  nJ.  (3) 

z  c  ^  x  y  y  x  ) 

This  field  gives  rise  to  a  diffusion  in  the  magnetic  field 


3B  .  2  3B 
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where  the  term  on  the  right  hand  side  is  obtained  from  Eq.  (2)  and  the 
factor  of  2  in  front  comes  from  adding  the  complex  conjugate.  Also,  we 
have  changed  the  notation  slightly  and  now  the  index  i  denotes  the  i*"*1 
rational  surface.  Thus  the  stabilized  resistive  interchange  mode 
gives  rise  to  magnetic  diffusion  on  each  rational  surface.  If  the 
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different  mode  widths  overlap,  there  will  be  magnetic  diffusion  over 
the  entire  plasma. 

In  addition  to  magnetic  diffusion,  the  fluctuations  also  give 
rise  to  energy  transport.  The  total  energy  equation  for  the  plasma  is 


Q  +  2 


P 


J  E 


z  z 


(5) 


where  p  is  the  pressure  and  Ez  is  given  by  Eq.  (3).  Note  that  because 
E  =  0,  there  is  no  E*J  term  on  the  right  hand  side  of  Eq.  (5).  As  is 
apparent  from  Eq.  (5),  the  presence  of  the  fluctuating  velocity  has 
two  effects.  First  it  gives  rise  to  an  anomalous  energy  flux  from 
the  —  p  term,  and  second  it  gives  rise  to  an  anomalous  energy 

exchange  with  the  magnetic  field  arising  from  the  (-Jz/c) (V^B^-V^B^  ) 
contribution  to  the  right  hand  side  of  Eq.  (5).  A  more  detailed 
discussion  of  the  energy  flux  will  be  published  elsewhere. 

To  get  a  rough  estimate  for  the  anomalous  magnetic  diffusion, 
let  us  assume  that  nonlinear  effects  limit  the  y  fluctuating  velocity 

-  X 

to  some  fraction  a  of  the  sound  speed  so  that  V  =  0.5  a  kH  r In 

o  M 

this  case,  the  value  of  the  magnetic  diffusion  coefficient  at  x  = 
is  given  roughly  by 


0  - 1  (&)  »2  *2  (s)  •  «> 

1/3 

It  is  interesting  to  note  that  in  this  case  D  q  as  one  might 

expect  for  resistive  interchange  modes.  Also  B  is  given  by 
V  x 

1/2  B  By  s&  1/2  BkH  a  BQ so  that  the  basic  phenomena  is  fluid 

convection  and  vortex  motion;  the  radial  fluctuating  magnetic  field 
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being,  by  contrast,  small.  For  Eta-Beta  II,  we  find  that  if  kr  ^  10, 
Eq.(6)  with  o<  1  gives  roughly  the  correct  confinement  time.  However, 
we  emphasize  that  since  l  ^  (the  ion  larmor  radius),  the  theory 
developed  here  does  not  directly  apply.  Nevertheless  it  would 
undoubtedly  be  interesting  to  try  to  measure  small  scale  velocity 
fluctuations  during  the  quiescent  phase  of  a  reverse  field  pinch. 

To  conclude,  we  have  shown  that  if  resistive  interchange  modes 
exist  in  a  nonlinearly  stabilized  state,  they  can  give  rise  to 
anomalous  magnetic  diffusion  which  is  localized  about  each  rational 
surface.  The  diffusion  coefficient  and  energy  transport  can  be 
calculated  easily  in  terms  of  the  velocity  fluctuation.  Since  simple 
theory  of  this  mode  shows  that  it  is  always  unstable  if  there  is  a 
pressure  gradient  opposite  in  sign  to  R,  it  seems  that  it  will  be  an 
Important  effect  in  both  RFP's  and  spheromaks. 

This  work  was  supported  by  the  Department  of  Energy.  The 
author  wishes  to  thank  Dr.  Joseph  DiMarco  for  providing  data  from 
Eta-Beta  II. 
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